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Abstract 

We analyze the Dalitz plot parameters of r\ — >3n decay in the framework of resummed chiral perturbation theory. This 
approach allows us to keep the uncertainties in the NNLO and higher orders under better control and estimate their 
influence. We cannot confirm the suspected discrepancy in the case of the charged decay parameter b, where even 
small uncertainties in higher orders could accommodate the difference. On the other hand, we find the experimental 
value of the neutral decay parameter a incompatible with an assumption of good convergence properties in the center 
of the Dalitz plot. We calculate nn rescattering bubble corrections up to three loops and show that these might explain 
the discrepancy, especially for a low value of the pseudoscalar decay constant in the chiral limit. However, that could 
indicate a failure of convergence of the chiral series in this channel already at low energies around 500MeV. 
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Chiral perturbation theory 0] 13 has a long history 
1 3 1 of trying to explain experimental data on the 77 —>3n 
decay. As is clear from the decay rate calculations, the 
theory converges really slowly for this decay channel. 
While the latest NNLO calculations J4) provide rea- 
sonable predictions for this quantity, experimental data 
are being gathered with increasing precision in order to 
make more detailed analysis of the Dalitz plot distribu- 
tion possible. The conventional Dalitz plot parameters 
are defined as (we generally follow notations from [4|): 

77 -> ttVtT : \A\ 2 = A 2 (l+ay + by 2 +dx 2 + ...) (1) 

77^371°: \A\ 2 =A 2 (l+az+...) (2) 

where x ~ u—t, y ~ sq-s, z ~ x 2 +y 2 and so is the Dalitz 
plot center s = 1 /3(M 2 +3M 2 ) . 

Comparison of the recent experimental information 
with the NNLO^PT results can be found in tables [TJ [2] 
As can be seen, the most obvious discrepancy appears 
in the charged decay parameter b and the neutral decay 
parameter a. 

Other approaches were developed in order to model 
the amplitude better, namely dispersive approaches lfT4l 
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a 


1998 Crystal Barrel (7) 
2001 Crystal Ball (8) 
2007 WASA at CELSIUS E) 
2009 WASA at COSY (TO) 
2009 Cr.Ball at MAMI-B (TT) 

2009 Cr.Ball at MAMI-C (12) 

2010 KLOE (TO) 


-0.052 ± 0.020 
-0.031 ±0.004 
-0.026 ± 0.014 
-0.027 ± 0.009 
-0.032 ± 0.003 
-0.0322 ± 0.0025 
-0.0301 ± 0.0050 


2007 NNLO ^PT @ 


+0.013 ± 0.032 



Table 1: Recent experimental and^PT results for the neutral channel. 





1998 2008 
Cr.Barrel(5) KLOE (6) 


2008 
NNLO ^PT (4) 


a 
b 
d 


-1.22 ±0.07 -1.090 ±0.020 
0.22 ±0.11 0.124 ±0.012 
0.06 ±0.04 0.057 ±0.017 


-1.271 ±0.075 
0.394 ±0.102 
0.055 ± 0.057 



Table 2: Recent experimental and^PT results for the charged channel. 
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Q21IT61Q2) and non-relativistic effective field theory ifTSl 
[l9l|20). These more or less abandon strict equivalence 
to ^PT and their success in reproducing a negative sign 
for a can serve as a motivation to ask what is the culprit 
of the failure of chiral perturbation theory to do so. 

Our aim is thus not to produce another alternative 
method, but to try to understand whether the theory, by 
which we mean ^PT as a low energy representation of 
QCD, really does have difficulties explaining the data 
and if so, try to identify the source of the problem. 

The starting point is the realization that the standard 
approach to ;fPT, as a usual treatment of perturbation 
series, implicitly assumes good convergence properties 
and hides the uncertainties associated with a possible 
violation of this assumption. Error bars are often not re- 
ported and it should be stressed that even the ones cited 
in the NNLO ^PT result above [4| are not systematic 
uncertainties inherent in the theory, but rather fitting er- 
rors tied to numerical procedure peculiar to the method 
used by the authors. 

There is a long standing suspicion that chiral per- 
turbation theory might posses a slow or irregular con- 
vergence in the case of the three quark flavor series 
ET1 l22l . the 77— >37r decay rate might serve as a prime 
example. An alternative approach, dubbed resummed 
^PT [23] [24), was developed in order to express these 
assumptions in terms of parameters and uncertainty 
bands. The procedure can be very shortly summarized 
in the following way: 

- standard ^PT Lagrangian and power counting 

- only expansions derived directly from the 
generating functional trusted 

- explicitly to NLO, higher orders collected 
in remainders 

- remainders retained, treated as sources of error 

- manipulations in non-perturbative algebraic way 

Our calculation closely follows the procedure out- 
lined in 11251 . What we present here is only a brief ex- 
cerpt, skipping all the details and concentrating only on 
the cases of the Dalitz plot parameters b and a. A more 
comprehensive work is in preparation ll32l . 

Within the formalism, we start by expressing the 
charged decay amplitude in terms of the 4-point Green 
functions Giju- We compute at first order in isospin 
braking, in this case the amplitude takes the form 

FlF n A( s ,t,u) = G 4 _ ra -e»G + _33 + e„G + _ 8 g + Ag, (3) 

where A^? is the direct higher order remainder to the 
complete 4-point Green function. The physical mixing 



angles to all chiral orders and first in isospin braking 
can be expressed in terms of quadratic mixing terms of 
the generating functional to NLO and related indirect 
remainders 

In this approximation the neutral decay channel ampli- 
tude can be related to the charged one as 

A(s, t, u) = A(s, t, u) + A(t, u, s) + A(u, s, t). (5) 

In accord with the method, strictly 0(p 2 ) parameters 
appear inside loops, while physical quantities in outer 
legs. Due to the leading order masses in loops such 
a strictly derived amplitude has an incorrect analytical 
structure, cuts and poles being in unphysical places. To 
account for this, the amplitude is carefully modified us- 
ing a NLO dispersive representation, the procedure is 
described in detail in [25 1. 

The next step is the treatment of the low energy con- 
stants (LEC's). The leading order ones are expressed in 
terms of convenient parameters 

Z = -4, X= , , , r=-r, R = -M '-, (6) 

F„ F£M£ m (m d -m u ) 

where m=(m„+m</)/2. The standard approach tacitly as- 
sumes values of X and Z close to one and r ~25, which 
means that the leading order terms should dominate the 
expansion. However, even the most recent standard 
^PT fit [26 1 indicates a much lower values of X and Z, 
thereby allowing for a possibility of a non-standard sce- 
nario of spontaneous chiral symmetry braking (SB^-S). 

At next-to-leading order, the LEC's L4-L8 are al- 
gebraically reparametrized in terms of pseudoscalat 
masses, decay constants and the free parameters X, Z 
and r using chiral expansions of two point Green func- 
tions, similarly to ll23l . Because expansions are for- 
mally not truncated, each generates an unknown higher 
order remainder. 

We don't have a similar procedure ready for Z4-L3 at 
this point, therefore we collect a set of standard ^PT fits 
[26 , 27 28 1 and by taking their mean and spread, while 
ignoring the much smaller reported error bars, we obtain 
an estimate of their influence. As will be shown in [32 1, 
the results depend on these constants only very weakly. 
The error bands given here include the estimated uncer- 
tainties in L1-L3. 

The 0(p b ) and higher order LEC's, notorious for their 
abundance, are collected in a relatively smaller num- 
ber of higher order remainders. We also fix the s-quark 
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Figure 1: Charged decay parameter b: top Z=l, bottom Z=0. 5 
horizontal lines: dashed - KLOE measurment |6|; dotted - S^PT |4] 
light band - statistical remainder estimate; dark band - result including 
nn rescattering, depending on scale /j=0.5^1GeV. 



Figure 2: Neutral decay parameter a: top Z=l, bottom Z=0.5 
horizontal lines: dashed - KLOE measurment [6]; dotted - S^PT |4| 
light band - statistical remainder estimate; dark band - result including 
TTTr rescattering, depending on scale //=0.5+lGeV. 



mass at r-25, motivated by lattice [29], as its value is 
anyway connected with X through the Kaplan-Manohar 
ambiguity [30 1 . We also investigated a low value r=15, 
but the results does not radically alter the presented con- 
clusions ll32l . At last, because at first order in isospin 
breaking the Dalitz plot parameters do not depend on 
R, we are left, besides the remainders, with two free 
parameters: X and Z. These effectively control the sce- 
nario of SByS in our results. 

Finally, the last step leading to numerical results is 
the estimate of the remainders. We have a direct re- 
mainder to the 4-point Green function and eight indirect 
ones - three related to both pseudoscalar masses and de- 
cay constants, two to mixing angles. We use two ap- 
proaches. The first one is based on general arguments 
about the convergence of the chiral series ll23l . which 
leads to 

A<? ~ ±0.1G, (7) 

where G here stands for any of our 2-point or 4-point 
Green functions, which generate the remainders. This 
is in principle an assumption. Hence we test the com- 
patibility of this assumption of a reasonably good chi- 
ral convergence of trusted quantities with experimen- 
tal data in a statistical sense. One peculiarity should 
be noted though - Dalitz plot parameters are deriva- 
tive quantities in terms of the Mandelstam variables and 
even if the direct remainder had a small absolute value 
around the center of the Dalitz plot, it could generate 
a large correction in derivatives. When making an ex- 
pansion of the direct remainder around the Dalitz plot 



center, such a large term would either grow to gener- 
ate a large direct remainder in some other region, possi- 
bly unphysical, but still where ^PT is assumed to con- 
verge well, or could be canceled by terms proportional 
to higher order derivatives, thus inducing some kind of 
fluctuation in the amplitude. As we feel such a behavior 
warrants explanation, we include the absence of unusu- 
ally large corrections in derivatives in our definition of 
good convergence properties and test for such a possi- 
bility as well. 

The results for the statistical remainder estimate are 
depicted in fig|T]and|2] Several conclusion can be made 
- in both cases the NNLO standard ^PT result lies in 
our uncertainty bands (the lighter ones), which is an im- 
portant consistency check. For the charged decay pa- 
rameter b, this is also true for the latest experimental 
measurement. This means we cannot confirm any dis- 
crepancy, as even a small correction compatible with the 
assumption reasonable chiral convergence could explain 
the differences between theory and experiment. On the 
other hand, the situation is quite different in the case of 
the neutral decay parameter a, where we can conclude, 
in accord with the NNLO S^PT result, that our defini- 
tion of good chiral convergence is not compatible with 
the data. It is also notable that this is true for any value 
of the free parameters X and Z and hence an alterna- 
tive scenario of SB^B, e.g. a small value of the quark 
condensate, is not the culprit here. 

The framework of resummed ^PT is well suited to 
include additional information about higher orders from 
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various sources l25l [3T1 . An independent estimate of 
the remainders can on one side be an important check 
of the validity of the statistical remainder estimate, as in 
the case of the parameter b, or could try to explain any 
deviances from this assumption, which is our aim for a. 

While other estimates being in preparation |32|, here 
we employ a specific higher order calculation, namely 
n-loop bubble contributions to final state nn rescatter- 
ing, which generally take the form 

n-i 

G£° = J] ^- k \s,Uu)£j' nn (s) l +(t,u channels), (8) 

k+l=\ 

where p„ and J^is) are the usual chiral logs and one 
loop scalar functions [2|, respectively. PVT \s,t,u) is 
an m-th order polynomial in the Mandelstam variables. 
We compute up to 0(p & ), that means 3-loop diagrams 
with LO vertices and 2-loop ones with one or two NLO 
counter terms. We stress that what we do is not a unita- 
rization procedure but a genuine ^-PT calculation. The 
motivation is that this is one of the suspects that could 
explain the discrepancy in a fl6l l20l [T7l . Terms with 
highest power in s are of the form ~N-jf^J(s,m JI ) n ~ l , 

where N is a numerical factor. In our case we obtain 
Ar(i-ioop) =1/2) Ar (2-ioop) = 4 /3) jv (3 - loop) =5/8, which im- 
plies that N is not a suppression factor with LO ver- 
tices. Thus at \^jj(s,m„) «1 convergence blows up, 

which a simple analysis can show is around ~600- 
700MeV at Z =1 or V* ~400-500MeV at Z=0.5, where 
the renormalization scale runs through /i=0.5-^lGeV. 

The result can be seen as the dark bands in figures 
[T] and [2] If they did not wander outside the light ones, 
it would indicate a confirmation that the nn rescattering 
bubble contributions agree with the statistical remainder 
estimate. This is indeed true for b, except a small area of 
the parameter space. Once again the case of a is differ- 
ent and we can see that kji rescattering could generate a 
negative sign, especially in a case of small value of the 
pseudoscalar decay constant in the chiral limit. 

Of course, this is very far from a complete calculation 
to <9(p 8 ), which also expresses itself in dependence on 
the renormalization scale p. But the fact that the scale 
dependence is quite benign could be also interpreted in 
the way that it is not unreasonable to consider the bubble 
nn rescattering separately. 

A further discussion can be made [32 1 and we will 
only summarize the results briefly - the source of the 
large negative 7T7r-rescattering contributions to a turns 
out to be where suspected, the terms leading in s. 
These generate a concavity in the amplitude, which is 
measured precisely by a through the second derivative. 



It is interesting to note that the 0(p & ) contributions to 
a are actually larger than the 0(p b ) ones, and a quick 
check of the anticipated form of the s-leading terms in 
even higher order contributions show that the next few 
orders can be expected to be large and possibly negative 
as well. These terms do not induce a large correction 
to the direct remainder in the center of the Dalitz plot, 
but the concavity is connected with a quick failure 
of the convergence of chiral series in the suspected 
energy region. That is unphysical for the case of the 
r] ->3tt decay, but could be an indication, if some other 
contribution do not cancel them, of a breakdown of 
the chiral expansion at quite low energies. This could 
be due to some higher energy structure present, for 
example a resonance. 
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